Construction of relatedness matrices in autopolyploid populations using low-depth high-throughput sequencing data.

KEY MESSAGE: An improved estimator of genomic relatedness using low-depth high-throughput sequencing data for autopolyploids is developed. Its outputs strongly correlate with SNP array-based estimates and are available in the package GUSrelate. High-throughput sequencing (HTS) methods have reduced sequencing costs and resources compared to array-based tools, facilitating the investigation of many non-model polyploid species. One important quantity that can be computed from HTS data is the genetic relatedness between all individuals in a population. However, HTS data are often messy, with multiple sources of errors (i.e. sequencing errors or missing parental alleles) which, if not accounted for, can lead to bias in genomic relatedness estimates. We derive a new estimator for constructing a genomic relationship matrix (GRM) from HTS data for autopolyploid species that accounts for errors associated with low sequencing depths, implemented in the R package GUSrelate. Simulations revealed that GUSrelate performed similarly to existing GRM methods at high depth but reduced bias in self-relatedness estimates when the sequencing depth was low. Using a panel consisting of 351 tetraploid potato genotypes, we found that GUSrelate produced GRMs from genotyping-by-sequencing (GBS) data that were highly correlated with a GRM computed from SNP array data, and less biased than existing methods when benchmarking against the array-based GRM estimates. GUSrelate provides researchers with a tool to reliably construct GRMs from low-depth HTS data.

Characterization of the Ruminal Microbiome of Water Buffaloes (Bubalus bubalis) Kept in Different Ecosystems in the Eastern Amazon.

Increasing the efficiency of rumen fermentation is one of the main ways to maximize the production of ruminants. It is therefore important to understand the ruminal microbiome, as well as environmental influences on that community. However, there are no studies that describe the ruminal microbiota in buffaloes in the Amazon. The objective of this study was to characterize the rumen microbiome of the water buffalo (Bubalus bubalis) in the eastern Amazon in the dry and rainy seasons in three grazing ecosystems: Baixo Amazonas (BA), Continente do Pará (CP), Ilha do Marajó (IM), and in a confinement system: Tomé-Açu (TA). Seventy-one crossbred male buffaloes (Murrah × Mediterranean) were used, aged between 24 and 36 months, with an average weight of 432 kg in the rainy season and 409 kg in the dry season, and fed on native or cultivated pastures. In the confinement system, the feed consisted of sorghum silage, soybean meal, wet sorghum premix, and commercial feed. Samples of the diet from each ecosystem were collected for bromatological analysis. The collections of ruminal content were carried out in slaughterhouses, with the rumen completely emptied and homogenized, the solid and liquid fractions separated, and the ruminal pH measured. DNA was extracted from the rumen samples, then sequenced using Restriction Enzyme Reduced Representation Sequencing. The taxonomic composition was largely similar between ecosystems. All 61 genera in the reference database were recognized, including members of the domains Bacteria and Archaea. The abundance of 23 bacterial genera differed significantly (p < 0.01) between the Tomé-Açu confinement and other ecosystems. Bacillus, Ruminococcus, and Bacteroides had lower abundance in samples from the Tomé-Açu system. Among the Archaea, the genus Methanomicrobium was less abundant in Tomé-Açu, while Methanosarcina was more abundant. There was a difference caused by all evaluated factors, but the diet (available or offered) was what most influenced the ruminal microbiota.

Large-scale analysis of sheep rumen metagenome profiles captured by reduced representation sequencing reveals individual profiles are influenced by the environment and genetics of the host.

BACKGROUND: Producing animal protein while reducing the animal's impact on the environment, e.g., through improved feed efficiency and lowered methane emissions, has gained interest in recent years. Genetic selection is one possible path to reduce the environmental impact of livestock production, but these traits are difficult and expensive to measure on many animals. The rumen microbiome may serve as a proxy for these traits due to its role in feed digestion. Restriction enzyme-reduced representation sequencing (RE-RRS) is a high-throughput and cost-effective approach to rumen metagenome profiling, but the systematic (e.g., sequencing) and biological factors influencing the resulting reference based (RB) and reference free (RF) profiles need to be explored before widespread industry adoption is possible. RESULTS: Metagenome profiles were generated by RE-RRS of 4,479 rumen samples collected from 1,708 sheep, and assigned to eight groups based on diet, age, time off feed, and country (New Zealand or Australia) at the time of sample collection. Systematic effects were found to have minimal influence on metagenome profiles. Diet was a major driver of differences between samples, followed by time off feed, then age of the sheep. The RF approach resulted in more reads being assigned per sample and afforded greater resolution when distinguishing between groups than the RB approach. Normalizing relative abundances within the sampling Cohort abolished structures related to age, diet, and time off feed, allowing a clear signal based on methane emissions to be elucidated. Genus-level abundances of rumen microbes showed low-to-moderate heritability and repeatability and were consistent between diets. CONCLUSIONS: Variation in rumen metagenomic profiles was influenced by diet, age, time off feed and genetics. Not accounting for environmental factors may limit the ability to associate the profile with traits of interest. However, these differences can be accounted for by adjusting for Cohort effects, revealing robust biological signals. The abundances of some genera were consistently heritable and repeatable across different environments, suggesting that metagenomic profiles could be used to predict an individual's future performance, or performance of its offspring, in a range of environments. These results highlight the potential of using rumen metagenomic profiles for selection purposes in a practical, agricultural setting.

Combining host and rumen metagenome profiling for selection in sheep: prediction of methane, feed efficiency, production, and health traits.

BACKGROUND: Rumen microbes break down complex dietary carbohydrates into energy sources for the host and are increasingly shown to be a key aspect of animal performance. Host genotypes can be combined with microbial DNA sequencing to predict performance traits or traits related to environmental impact, such as enteric methane emissions. Metagenome profiles were generated from 3139 rumen samples, collected from 1200 dual purpose ewes, using restriction enzyme-reduced representation sequencing (RE-RRS). Phenotypes were available for methane (CH4) and carbon dioxide (CO2) emissions, the ratio of CH4 to CH4 plus CO2 (CH4Ratio), feed efficiency (residual feed intake: RFI), liveweight at the time of methane collection (LW), liveweight at 8 months (LW8), fleece weight at 12 months (FW12) and parasite resistance measured by faecal egg count (FEC1). We estimated the proportion of phenotypic variance explained by host genetics and the rumen microbiome, as well as prediction accuracies for each of these traits. RESULTS: Incorporating metagenome profiles increased the variance explained and prediction accuracy compared to fitting only genomics for all traits except for CO2 emissions when animals were on a grass diet. Combining the metagenome profile with host genotype from lambs explained more than 70% of the variation in methane emissions and residual feed intake. Predictions were generally more accurate when incorporating metagenome profiles compared to genetics alone, even when considering profiles collected at different ages (lamb vs adult), or on different feeds (grass vs lucerne pellet). A reference-free approach to metagenome profiling performed better than metagenome profiles that were restricted to capturing genera from a reference database. We hypothesise that our reference-free approach is likely to outperform other reference-based approaches such as 16S rRNA gene sequencing for use in prediction of individual animal performance. CONCLUSIONS: This paper shows the potential of using RE-RRS as a low-cost, high-throughput approach for generating metagenome profiles on thousands of animals for improved prediction of economically and environmentally important traits. A reference-free approach using a microbial relationship matrix from log10 proportions of each tag normalized within cohort (i.e., the group of animals sampled at the same time) is recommended for future predictions using RE-RRS metagenome profiles.

Genetic parameter analysis of bareness and tail traits in New Zealand sheep.

Physical traits that improve welfare and disease outcomes for sheep are becoming increasingly important due to both increased climate challenges and societal expectations. Such traits include tail length, the amount of skin (vs. wool) on the underside of the tail, and the area of no-wool (hair) on the belly and breech areas (surrounding the anus) of the animal. An industry dataset consisting of records from individual stud breeders and industry progeny tests was available to estimate the genetic parameters associated with these traits and to investigate the potential for within-breed genetic selection. The heritability estimate for tail length was 0.68 ± 0.01 when breed was not fitted, and 0.63 ± 0.01 when breed was fitted. Similar trends were observed for breech and belly bareness which had heritability estimates around 0.50 (± 0.01). The estimates for these bareness traits are both higher than previous reports from animals of the same age. There was, however, between breed variation in the starting point for these traits, with some breeds having significantly longer tails and a wooly breech and belly, and limited variability. Overall, the results of this study show that flocks exhibiting some variation will be able to make rapid genetic progress in selecting for bareness and tail length traits, and therefore have the potential to make progress towards a sheep that is easier to look after and suffers fewer welfare insults. For those breeds that showed limited within-breed variation, outcrossing may be required to introduce genotypes that exhibit shorter tail length and bareness of belly and breech to increase the rate of genetic gain. Whatever approach is taken by the industry, these results support that genetic improvement can be used to breed "ethically improved sheep".

Physical traits that improve welfare and disease outcomes for sheep are becoming increasingly important as a result of both increased climate challenges and societal expectations. Such physical traits include genetically short tails and bare (hair vs. wool) in the breech and belly region. This study demonstrates that these traits are highly heritable, indicating that a lot of the variation observed between individuals is due to genetic variation. Overall, the results of this study show that flocks exhibiting some variation will be able to make rapid genetic progress in selecting for bareness and tail length traits, and therefore make progress towards a sheep that is easier to look after and suffers fewer welfare insults. Some breeds, however, showed very little within-breed variation, and outcrossing may be required to introduce genotypes that exhibit the traits. Whatever approach is taken by the industry, these results support that genetic improvement can be used to breed "ethically improved sheep".

Efficiency of genotyping by sequencing in inferring genomic relatedness and molecular insights into fat tail selection in Tunisian sheep.

In developing countries, the use of simple and cost-efficient molecular technology is crucial for genetic characterization of local animal resources and better development of conservation strategies. The genotyping by sequencing (GBS) technique, also called restriction enzyme- reduced representational sequencing, is an efficient, cost-effective method for simultaneous discovery and genotyping of many markers. In the present study, we applied a two-enzyme GBS protocol (PstI/MspI) to discover and genotype SNP markers among 197 Tunisian sheep samples. A total of 100 333 bi-allelic SNPs were discovered and genotyped with an SNP call rate of 0.69 and mean sample depth 3.33. The genomic relatedness between 183 samples grouped the samples perfectly to their populations and pointed out a high genetic relatedness of inbred subpopulation reflecting the current adopted reproductive strategies. The genome-wide association study contrasting fat vs. thin-tailed breeds detected 41 significant variants including a peak positioned on OAR20. We identified FOXC1, GMDS, VEGFA, OXCT1, VRTN and BMP2 as the most promising for sheep tail-type trait. The GBS data have been useful to assess the population structure and improve our understanding of the genomic architecture of distinctive characteristics shaped by selection pressure in local sheep breeds. This study successfully investigates a cost-efficient method to discover genotypes, assign populations and understand insights into sheep adaptation to arid area. GBS could be of potential utility in livestock species in developing/emerging countries.

Impact of breeding for reduced methane emissions in New Zealand sheep on maternal and health traits.

Enteric methane emissions from ruminants account for ∼35% of New Zealand's greenhouse gas emissions. This poses a significant threat to the pastoral sector. Breeding has been shown to successfully lower methane emissions, and genomic prediction for lowered methane emissions has been introduced at the national level. The long-term genetic impacts of including low methane in ruminant breeding programs, however, are unknown. The success of the New Zealand sheep industry is currently heavily reliant on the prolificacy, fecundity and survival of adult ewes. The objective of this study was to determine genetic and phenotypic correlations between adult maternal ewe traits (live weight, body condition score, number of lambs born, litter survival to weaning, pregnancy scanning and fleece weight), faecal and Nematodirus egg counts and measures of methane in respiration chambers. More than 9,000 records for methane from over 2,200 sheep measured in respiration chambers were collected over 10 years. Sheep were fed on a restricted diet calculated as approximately twice the maintenance. Methane measures were converted to absolute daily emissions of methane measured in g per day (CH4/day). Two measures of methane yield were recorded: the ratio of CH4 to dry matter intake (g CH4/kg DMI; CH4/DMI) and the ratio of CH4 to total gas emissions (CH4/(CH4 + CO2)). Ewes were maintained in the flocks for at least two parities. Non-methane trait data from over 8,000 female relatives were collated to estimate genetic correlations. Results suggest that breeding for low CH4/DMI is unlikely to negatively affect faecal egg counts, adult ewe fertility and litter survival traits, with no evidence for significant genetic correlations. Fleece weight was unfavourably (favourably) correlated with CH4/DMI (rg = -0.21 ± 0.09). Live weight (rg = 0.3 ± 0.1) and body condition score (rg = 0.2 ± 0.1) were positively correlated with methane yield. Comparing the two estimates of methane yield, CH4/DMI had lower heritability and repeatability. However, correlations of both measures with adult ewe traits were similar. This suggests that breeding is a suitable mitigation strategy for lowering methane yield, but wool, live weight and fat deposition traits may be affected over time and should be monitored.

Genetic parameters for residual feed intake, methane emissions, and body composition in New Zealand maternal sheep.

There is simultaneous interest in improving the feed efficiency of ruminant livestock and reducing methane (CH4) emissions. The relationship (genetic and phenotypic) between feed efficiency (characterized as residual feed intake: RFI) and greenhouse gases [methane (CH4) and carbon dioxide (CO2)] traits in New Zealand (NZ) maternal sheep has not previously been investigated, nor has their relationship with detailed estimates of body composition. To investigate these relationships in NZ maternal sheep, a feed intake facility was established at AgResearch Invermay, Mosgiel, NZ in 2015, comprising automated feeders that record individual feeding events. Individual measures of feed intake, feeding behavior (length and duration of eating events), and gas emissions (estimated using portable accumulation chambers) were generated on 986 growing maternal ewe lambs sourced from three pedigree recorded flocks registered in the Sheep Improvement Limited database (www.sil.co.nz). Additional data were generated from a subset of 591 animals for body composition (estimated using ultrasound and computed tomography scanning). The heritability estimates for RFI, CH4, and CH4/(CH4+CO2) were 0.42 ± 0.09, 0.32 ± 0.08, and 0.29 ± 0.06, respectively. The heritability estimates for the body composition traits were high for carcass lean and fat traits; for example, the heritability for visceral fat (adjusted for body weight) was 0.93 ± 0.19. The relationship between RFI and CH4 emissions was complex, and although less feed eaten will lead to a lowered absolute amount of CH4 emitted, there was a negative phenotypic and genetic correlation between RFI and CH4/(CH4+CO2) of -0.13 ± 0.03 and -0.41 ± 0.15, respectively. There were also genetic correlations, that were different from zero, between both RFI and CH4 traits with body composition including a negative correlation between the proportion of visceral fat in the body and RFI (-0.52 ± 0.16) and a positive correlation between the proportion of lean in the body and CH4 (0.54 ± 0.12). Together the results provide the first accurate estimates of the genetic correlations between RFI, CH4 emissions, and the body composition (lean and fat) in sheep. These correlations will need to be accounted for in genetic improvement programs.

Low-cost sample preservation methods for high-throughput processing of rumen microbiomes.

BACKGROUND: The use of rumen microbial community (RMC) profiles to predict methane emissions has driven interest in ruminal DNA preservation and extraction protocols that can be processed cheaply while also maintaining or improving DNA quality for RMC profiling. Our standard approach for preserving rumen samples, as defined in the Global Rumen Census (GRC), requires time-consuming pre-processing steps of freeze drying and grinding prior to international transportation and DNA extraction. This impedes researchers unable to access sufficient funding or infrastructure. To circumvent these pre-processing steps, we investigated three methods of preserving rumen samples for subsequent DNA extraction, based on existing lysis buffers Tris-NaCl-EDTA-SDS (TNx2) and guanidine hydrochloride (GHx2), or 100% ethanol. RESULTS: Rumen samples were collected via stomach intubation from 151 sheep at two time-points 2 weeks apart. Each sample was separated into four subsamples and preserved using the three preservation methods and the GRC method (n = 4 × 302). DNA was extracted and sequenced using Restriction Enzyme-Reduced Representation Sequencing to generate RMC profiles. Differences in DNA yield, quality and integrity, and sequencing metrics were observed across the methods (p < 0.0001). Ethanol exhibited poorer quality DNA (A260/A230 < 2) and more failed samples compared to the other methods. Samples preserved using the GRC method had smaller relative abundances in gram-negative genera Anaerovibrio, Bacteroides, Prevotella, Selenomonas, and Succiniclasticum, but larger relative abundances in the majority of 56 additional genera compared to TNx2 and GHx2. However, log10 relative abundances across all genera and time-points for TNx2 and GHx2 were on average consistent (R2 > 0.99) but slightly more variable compared to the GRC method. Relative abundances were moderately to highly correlated (0.68 ± 0.13) between methods for samples collected within a time-point, which was greater than the average correlation (0.17 ± 0.11) between time-points within a preservation method. CONCLUSIONS: The two modified lysis buffers solutions (TNx2 and GHx2) proposed in this study were shown to be viable alternatives to the GRC method for RMC profiling in sheep. Use of these preservative solutions reduces cost and improves throughput associated with processing and sequencing ruminal samples. This development could significantly advance implementation of RMC profiles as a tool for breeding ruminant livestock.

Hitchhiking Mapping of Candidate Regions Associated with Fat Deposition in Iranian Thin and Fat Tail Sheep Breeds Suggests New Insights into Molecular Aspects of Fat Tail Selection.

The fat tail is a phenotype that divides indigenous Iranian sheep genetic resources into two major groups. The objective of the present study is to refine the map location of candidate regions associated with fat deposition, obtained via two separate whole genome scans contrasting thin and fat tail breeds, and to determine the nature of the selection occurring in these regions using a hitchhiking approach. Zel (thin tail) and Lori-Bakhtiari (fat tail) breed samples that had previously been run on the Illumina Ovine 50 k BeadChip, were genotyped with a denser set of SNPs in the three candidate regions using a Sequenom Mass ARRAY platform. Statistical tests were then performed using different and complementary methods based on either site frequency (FST and Median homozygosity) or haplotype (iHS and XP-EHH). The results from candidate regions on chromosome 5 and X revealed clear evidence of selection with the derived haplotypes that was consistent with selection to near fixation for the haplotypes affecting fat tail size in the fat tail breed. An analysis of the candidate region on chromosome 7 indicated that selection differentiated the beneficial alleles between breeds and homozygosity has increased in the thin tail breed which also had the ancestral haplotype. These results enabled us to confirm the signature of selection in these regions and refine the critical intervals from 113 kb, 201 kb, and 2831 kb to 28 kb, 142 kb, and 1006 kb on chromosome 5, 7, and X respectively. These regions contain several genes associated with fat metabolism or developmental processes consisting of TCF7 and PPP2CA (OAR5), PTGDR and NID2 (OAR7), AR, EBP, CACNA1F, HSD17B10,SLC35A2, BMP15, WDR13, and RBM3 (OAR X), and each of which could potentially be the actual target of selection. The study of core haplotypes alleles in our regions of interest also supported the hypothesis that the first domesticated sheep were thin tailed, and that fat tail animals were developed later. Overall, our results provide a comprehensive assessment of how and where selection has affected the patterns of variation in candidate regions associated with fat deposition in thin and fat tail sheep breeds.

Genetic variation in skin traits in New Zealand lambs.

BACKGROUND: This study explored the genetic variability in the New Zealand sheep population for economically important skin traits. Skins were collected at slaughter from two progeny test flocks, resulting in 725 skins evaluated for grain strain, flatness, crust leather strength and overall suitability for shoe leather. DNA profiles collected from skins post-slaughter were matched to individual animals using previously collected high-density genotypes. RESULTS: Considerable phenotypic variation for skin traits was observed, with around 40% of the skins being identified as suitable for high-value shoe leather production. Several key traits associated with leather production, including flatness, tear strength, grain strength and grain strain were found to be moderate to highly heritable (h2 = 0.28-0.82). There were no major significant genome-wide association study (GWAS) peaks associated with many of the traits examined, however, one single-nucleotide polymorphism (SNP) reached significance for the flatness of the skin over the hindquarters. CONCLUSION: This research confirms that suitable lamb skins can be bred for use as high-value shoe leather. While moderately to highly heritable, skin traits in New Zealand lambs appear to be polygenic with no genes of major effect underlaying the traits of interest. Given the complex nature of these traits, the identification and selection of animals with higher-value skins may be enabled by geomic selection. © 2022 The Authors. Journal of The Science of Food and Agriculture published by John Wiley & Sons Ltd on behalf of Society of Chemical Industry.

Genomic Tools for the Identification of Loci Associated with Facial Eczema in New Zealand Sheep.

Facial eczema (FE) is a significant metabolic disease that affects New Zealand ruminants. Ingestion of the mycotoxin sporidesmin leads to liver and bile duct damage, which can result in photosensitisation, reduced productivity and death. Strategies used to manage the incidence and severity of the disease include breeding. In sheep, there is considerable genetic variation in the response to FE. A commercial testing program is available for ram breeders who aim to increase tolerance, determined by the concentration of the serum enzyme, gamma-glutamyltransferase 21 days after a measured sporidesmin challenge (GGT21). Genome-wide association studies were carried out to determine regions of the genome associated with GGT21. Two regions on chromosomes 15 and 24 are reported, which explain 5% and 1% of the phenotypic variance in the response to FE, respectively. The region on chromosome 15 contains the ß-globin locus. Of the significant SNPs in the region, one is a missense variant within the haemoglobin subunit ß (HBB) gene. Mass spectrometry of haemoglobin from animals with differing genotypes at this locus indicated that genotypes are associated with different forms of adult ß-globin. Haemoglobin haplotypes have previously been associated with variation in several health-related traits in sheep and warrant further investigation regarding their role in tolerance to FE in sheep. We show a strategic approach to the identification of regions of importance for commercial breeding programs with a combination of discovery, statistical and biological validation. This study highlights the power of using increased density genotyping for the identification of influential genomic regions, combined with subsequent inclusion on lower density genotyping platforms.

Genomic signatures of inbreeding in a critically endangered parrot, the kakapo.

Events of inbreeding are inevitable in critically endangered species. Reduced population sizes and unique life-history traits can increase the severity of inbreeding, leading to declines in fitness and increased risk of extinction. Here, we investigate levels of inbreeding in a critically endangered flightless parrot, the kakapo (Strigops habroptilus), wherein a highly inbred island population and one individual from the mainland of New Zealand founded the entire extant population. Genotyping-by-sequencing (GBS), and a genotype calling approach using a chromosome-level genome assembly, identified a filtered set of 12,241 single-nucleotide polymorphisms (SNPs) among 161 kakapo, which together encompass the total genetic potential of the extant population. Multiple molecular-based estimates of inbreeding were compared, including genome-wide estimates of heterozygosity (FH), the diagonal elements of a genomic-relatedness matrix (FGRM), and runs of homozygosity (RoH, FRoH). In addition, we compared levels of inbreeding in chicks from a recent breeding season to examine if inbreeding is associated with offspring survival. The density of SNPs generated with GBS was sufficient to identify chromosomes that were largely homozygous with RoH distributed in similar patterns to other inbred species. Measures of inbreeding were largely correlated and differed significantly between descendants of the two founding populations. However, neither inbreeding nor ancestry was found to be associated with reduced survivorship in chicks, owing to unexpected mortality in chicks exhibiting low levels of inbreeding. Our study highlights important considerations for estimating inbreeding in critically endangered species, such as the impacts of small population sizes and admixture between diverse lineages.

Development of Epigenetic Clocks for Key Ruminant Species.

Robust biomarkers of chronological age have been developed in humans and model mammalian species such as rats and mice using DNA methylation data. The concept of these so-called "epigenetic clocks" has emerged from a large body of literature describing the relationship between genome-wide methylation levels and age. Epigenetic clocks exploit this phenomenon and use small panels of differentially methylated cytosine (CpG) sites to make robust predictions of chronological age, independent of tissue type. Here, we present highly accurate livestock epigenetic clocks for which we have used the custom mammalian methylation array "HorvathMammalMethyl40" to construct the first epigenetic clock for domesticated goat (Capra hircus), cattle (Bos taurus), Red (Cervus elaphus) and Wapiti deer (Cervus canadensis) and composite-breed sheep (Ovis aries). Additionally, we have constructed a 'farm animal clock' for all species included in the study, which will allow for robust predictions to be extended to various breeds/strains. The farm animal clock shows similarly high accuracy to the individual species' clocks (r > 0.97), utilizing only 217 CpG sites to estimate age (relative to the maximum lifespan of the species) with a single mathematical model. We hypothesise that the applications of this livestock clock could extend well beyond the scope of chronological age estimates. Many independent studies have demonstrated that a deviation between true age and clock derived molecular age is indicative of past and/or present health (including stress) status. There is, therefore, untapped potential to utilize livestock clocks in breeding programs as a predictor for age-related production traits.

Application of Low Coverage Genotyping by Sequencing in Selectively Bred Arctic Charr (Salvelinus alpinus).

Arctic charr (Salvelinus alpinus) is a species of high economic value for the aquaculture industry, and of high ecological value due to its Holarctic distribution in both marine and freshwater environments. Novel genome sequencing approaches enable the study of population and quantitative genetic parameters even on species with limited or no prior genomic resources. Low coverage genotyping by sequencing (GBS) was applied in a selected strain of Arctic charr in Sweden originating from a landlocked freshwater population. For the needs of the current study, animals from year classes 2013 (171 animals, parental population) and 2017 (759 animals; 13 full sib families) were used as a template for identifying genome wide single nucleotide polymorphisms (SNPs). GBS libraries were constructed using the PstI and MspI restriction enzymes. Approximately 14.5K SNPs passed quality control and were used for estimating a genomic relationship matrix. Thereafter a wide range of analyses were conducted in order to gain insights regarding genetic diversity and investigate the efficiency of the genomic information for parentage assignment and breeding value estimation. Heterozygosity estimates for both year classes suggested a slight excess of heterozygotes. Furthermore, FST estimates among the families of year class 2017 ranged between 0.009 - 0.066. Principal components analysis (PCA) and discriminant analysis of principal components (DAPC) were applied aiming to identify the existence of genetic clusters among the studied population. Results obtained were in accordance with pedigree records allowing the identification of individual families. Additionally, DNA parentage verification was performed, with results in accordance with the pedigree records with the exception of a putative dam where full sib genotypes suggested a potential recording error. Breeding value estimation for juvenile growth through the usage of the estimated genomic relationship matrix clearly outperformed the pedigree equivalent in terms of prediction accuracy (0.51 opposed to 0.31). Overall, low coverage GBS has proven to be a cost-effective genotyping platform that is expected to boost the selection efficiency of the Arctic charr breeding program.

Exclusion and Genomic Relatedness Methods for Assignment of Parentage Using Genotyping-by-Sequencing Data.

Genotypes are often used to assign parentage in agricultural and ecological settings. Sequencing can be used to obtain genotypes but does not provide unambiguous genotype calls, especially when sequencing depth is low in order to reduce costs. In that case, standard parentage analysis methods no longer apply. A strategy for using low-depth sequencing data for parentage assignment is developed here. It entails the use of relatedness estimates along with a metric termed excess mismatch rate which, for parent-offspring pairs or trios, is the difference between the observed mismatch rate and the rate expected under a model of inheritance and allele reads without error. When more than one putative parent has similar statistics, bootstrapping can provide a measure of the relatedness similarity. Putative parent-offspring trios can be further checked for consistency by comparing the offspring's estimated inbreeding to half the parent relatedness. Suitable thresholds are required for each metric. These methods were applied to a deer breeding operation consisting of two herds of different breeds. Relatedness estimates were more in line with expectation when the herds were analyzed separately than when combined, although this did not alter which parents were the best matches with each offspring. Parentage results were largely consistent with those based on a microsatellite parentage panel with three discordant parent assignments out of 1561. Two models are investigated to allow the parentage metrics to be calculated with non-random selection of alleles. The tools and strategies given here allow parentage to be assigned from low-depth sequencing data.

Genetic parameters of plasma and ruminal volatile fatty acids in sheep fed alfalfa pellets and genetic correlations with enteric methane emissions1.

Animal-to-animal variation in methane (CH4) emissions determined in respiration chambers has a genetic basis, but rapid phenotyping methods that can be applied on-farm are required to enable increased genetic progress by the farming industry. Fermentation of carbohydrates in the rumen results in the formation of VFA with hydrogen (H2) as a byproduct that is used for CH4 formation. Generally, fermentation pathways leading to acetate are associated with the most H2 production, less H2 formation is associated with butyrate production, and propionate and valerate production are associated with reduced H2 production. Therefore, VFA may constitute a potential correlated proxy for CH4 emissions to enable high-throughput animal screening. The objective of the present study was to determine the genetic parameters for ruminal and plasma VFA concentrations in sheep fed alfalfa (Medicago sativa L.) pellets and their genetic (rg) and phenotypic (rp) correlations with CH4 emissions. Measurements of CH4 emissions in respiration chambers and ruminal (stomach tubing 18 h from last meal) and blood plasma (3 h post-feeding) VFA concentrations were made on 1,538 lambs from 5 birth years (2007 and 2009 to 2012) aged between 5 and 10 mo, while the animals were fed alfalfa pellets at 2.0 times maintenance requirements in 2 equal size meals (0900 and 1500 h). These measurements were repeated twice (rounds) 14 d apart. Mean (± SD) CH4 production was 24.4 ± 3.08 g/d, and the mean CH4 yield was 15.8 ± 1.51 g/kg DMI. Mean concentration of total ruminal VFA was 52.2 mM, with concentrations of acetate, propionate and butyrate of 35.97, 8.83, and 4.02 mM, respectively. Ruminal total VFA concentration had heritability (h2) and repeatability estimates (± SE) of 0.24 ± 0.05 and 0.35 ± 0.03, respectively, and similar estimates were found for acetate, propionate, and butyrate. Blood plasma concentrations of VFA had much lower estimates of h2 and repeatability than ruminal VFA. Genetic correlations with CH4 yield were greatest for total concentrations of ruminal VFA and acetate, with 0.54 ± 0.12 and 0.56 ± 0.12, respectively, which were much greater than their corresponding rp. The rp and rg of ruminal VFA proportions and blood VFAs with CH4 emissions were in general lower than for ruminal VFA concentrations. However, minor ruminal VFA proportions had also moderate rg with CH4 yield. Pre-feeding concentrations of total VFA and acetate were the strongest correlated proxies to select sheep that are genetically low CH4 emitters.

Developing Successful Breeding Programs for New Zealand Aquaculture: A Perspective on Progress and Future Genomic Opportunities.

Over the past 40 years New Zealand (NZ) aquaculture has grown into a significant primary industry. Tonnage is small on a global scale, but the industry has built an international reputation for the supply of high quality seafood to many overseas markets. Since the early 1990s the industry has recognized the potential gains from selective breeding and the challenge has been to develop programs that can overcome biological obstacles (such as larval rearing and mortality) and operate cost-effectively on a relatively small scale while still providing significant gains in multiple traits of economic value. This paper provides an overview of the current status, and a perspective on genomic technology implementation, for the family based genetic improvement programs established for the two main species farmed in NZ: Chinook (king) salmon (Oncorhynchus tshawytscha) and GreenshellTM mussel (Perna canaliculus). These programs have provided significant benefit to the industry in which we are now developing genomic resources based on genotyping-by-sequencing to complement the breeding programs, enable evaluation of the genetic diversity and identify the potential benefits of genomic selection. This represents an opportunity to increase genetic gain and more effectively utilize the potential for within family selection.

Reduced representation sequencing detects only subtle regional structure in a heavily exploited and rapidly recolonizing marine mammal species.

Next-generation reduced representation sequencing (RRS) approaches show great potential for resolving the structure of wild populations. However, the population structure of species that have shown rapid demographic recovery following severe population bottlenecks may still prove difficult to resolve due to high gene flow between subpopulations. Here, we tested the effectiveness of the RRS method Genotyping-By-Sequencing (GBS) for describing the population structure of the New Zealand fur seal (NZFS, Arctocephalus forsteri), a species that was heavily exploited by the 19th century commercial sealing industry and has since rapidly recolonized most of its former range from a few isolated colonies. Using 26,026 neutral single nucleotide polymorphisms (SNPs), we assessed genetic variation within and between NZFS colonies. We identified low levels of population differentiation across the species range (<1% of variation explained by regional differences) suggesting a state of near panmixia. Nonetheless, we observed subtle population substructure between West Coast and Southern East Coast colonies and a weak, but significant (p = 0.01), isolation-by-distance pattern among the eight colonies studied. Furthermore, our demographic reconstructions supported severe bottlenecks with potential 10-fold and 250-fold declines in response to Polynesian and European hunting, respectively. Finally, we were able to assign individuals treated as unknowns to their regions of origin with high confidence (96%) using our SNP data. Our results indicate that while it may be difficult to detect population structure in species that have experienced rapid recovery, next-generation markers and methods are powerful tools for resolving fine-scale structure and informing conservation and management efforts.

Genetic parameters of methane emissions determined using portable accumulation chambers in lambs and ewes grazing pasture and genetic correlations with emissions determined in respiration chambers.

Methane (CH4) emission traits were previously found to be heritable and repeatable in sheep fed alfalfa pellets in respiration chambers (RC). More rapid screening methods are, however, required to increase genetic progress and to provide a cost-effective method to the farming industry for maintaining the generation of breeding values in the future. The objective of the current study was to determine CH4 and carbon dioxide (CO2) emissions using several 1-h portable accumulation chamber (PAC) measurements from lambs and again as ewes while grazing ryegrass-based pasture. Many animals with PAC measurements were also measured in RC while fed alfalfa pellets at 2.0 × maintenance metabolizable energy requirements (MEm). Heritability estimates from mixed models for CH4 and CO2 production (g/d) were 0.19 and 0.16, respectively, when measured using PAC with lambs; 0.20 and 0.27, respectively, when measured using PAC with ewes; and 0.23 and 0.34, respectively, when measured using RC with lambs. For measured gas traits, repeatabilities of measurements collected 14 d apart ranged from 0.33 to 0.55 for PAC (combined lambs and ewes) and were greater at 0.65 to 0.76 for the same traits measured using RC. Genetic correlations (rg) between PAC in lambs and ewes were 0.99 for CH4, 0.93 for CH4 + CO2, and 0.85 for CH4/(CH4 + CO2), suggesting that CH4 emissions in lambs and ewes are the same trait. Genetic correlations between PAC and RC measurements were lower, at 0.62 to 0.67 for CH4 and 0.41 to 0.42 for CH4 + CO2, likely reflecting different environmental conditions associated with the protocols used with the 2 measurement methods. The CH4/(CH4 + CO2) ratio was the most similar genetic trait measured using PAC (both lambs and ewes, 63% and 66% selection efficiency, respectively) compared with CH4 yield (g/kg DMI) measured using RC. These results suggest that PAC measurements have considerable value as a rapid low-cost method to estimate breeding values for CH4 emissions in sheep.